A New Limit for the Coupling Constant of the Electron-Nucleus Scalar-Pseudoscalar 
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We report the results of our calculations of the atomic electric dipole moments of cesium and thallium 
arising from the electron-nucleus scalar-pseudoscalar interaction. The calculations are based on the 
all order relativistic coupled-cluster theory. Electron correlation effects, particularly in the case of 
thallium are of crucial importance. We obtain a new limit for the scalar-pseudoscalar interaction 
by combining the result of our thallium calculation and the measured value of the electric dipole 
moment of that atom. 



It is now widely recognized that atomic electric dipole 
moments (EDMs) arising from the violations of parity 
(P) and time-reversal (T) symmetries can provide im- 
portant information about new physics beyond the Stan- 
dard Model (SM) p. T violation implies CP violation 
via the CPT theorem . The dominant sources of the 
EDM of a paramagnetic atom are the EDM of an elec- 
tron and the scalar-pseudoscalar (S-PS) interaction be- 
tween the electron and the nucleus which violates P as 
well as T symmetries Q • While the former has been dealt 
with extensively Q , the latter has received relatively lit- 
tle attention. Barr has shown that in two Higgs doublet 
models 0, the contribution of the S-PS interaction to 
an atomic EDM can exceed that due to the EDM of an 
electron for certain values of the model parameters. The 
S-PS electron-nucleus interaction which originates from 
the S-PS electron-quark interactions could therefore shed 
light on new physics beyond the Standard Model. 

In this Letter we calculate the ratios of the atomic 
EDMs to the S-PS coupling constants for cesium (Cs) and 
thallium (Tl) using relativistic coupled-cluster (RCC) 
theory. Combining the result of our Tl calculation with 
experimental data we obtain a new limit for the S-PS cou- 
pling constant [Cs). Following the EDM measurement 
of atomic Cs by Murthy et al. 0, a series of experi- 
ments were carried out by Commins and co-workers on 
atomic Tl 0, H, • Currently the most accurate data on 
the measurement of the EDM of any paramagnetic atom 
comes from the experiment on Tl by Regan et al. 0- 



New experiments on the EDM of atomic Cs using the 
techniques of laser cooling and t rapping ar e in progress 
in three different laboratories pXllllLll3 . ll3 |. They could 
improve the existing limit of the atomic EDM by about 
two orders of magnitude. 

The P and T violating S-PS electron-nucleon interac- 
tion Hamiltonian (H|"q^) is given by 



S-PS 
EDM 



= ^CsAY,^(iellpN{Te). (0.1) 



In the above expression G_f is the Fermi constant, Cs is 
the dimensionless S-PS constant which is given by Cs — 
{ZCs,p + NCs,n)/A, where Z, N and A are the atomic, 
neutron and atomic mass numbers respectively. It is a 
weighted average of the electron-neutron and electron- 
proton coupling constants, Cs,n and Cs,p with C's.n ~ 
Cs,p- PNi^e) is the nuclear density function and 7^(= 
j7'^7^7^7^), which is a pseudo-scalar, is the product of 
the four Dirac matrices. 

Hedm responsible for mixing states of opposite pari- 
ties but with the same angular momentum. The strength 
of this interaction Hamiltonian is sufhciently weak for it 
to be considered as a first-order perturbation. It is there- 
fore possible to write the atomic wave function as 



(0.2) 



In RCC, the ground state |\E'[°h for a single valence (v) 
open-shell system is given hj 0, 



(0.3) 



tE-mail: bijaya@iiap.res.in, B.K.Sahoo@gsi.de 



where we define |<I>i,) — aJ^|<i>o), with |<I>o) as the Dirac- 
Fock (DF) state for closed-shell system. 
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In the singles and doubles approximation we have 



c(0) _ c(0) 



.(0) 



In the present work we obtain the first order wave func- 
tion as the solution of the equation 



-(0) 
'2d 



(0.4) 



(H(0)-i?("))|*«) = (i?W 



(0.7) 



where t|°^ and T^"^ are the single and double particle- where vanishes for first order correction. In this 
hole excitation operators for core electrons (holes) and 



and are the single and double excitation op- 
erators with the valence electron respectively. The am- 
plitudes corresponding to these operators can be deter- 
mined by solving the RCC singles and doubles equations 
(iGf. A subset of important triple excitations have been 
considered in the determination of the open shell ampli- 
tudes and which is described in 

Using eqn. (0.2) the EDM of an atom of the ground 
state is given by 



Da 



(^i°^|D|^i'V(^l^^|D|^l,°^) 
(*l°)|v|/(°') 



(0.5) 



where higher order perturbed terms are neglected and D 
is the electric dipole (El) operator. Using the explicit 
expression for the first order perturbed wave function 
(l^-l^^)), we get 



Da 



vI'f|D|vE-^"^)(*riH|D^^|*i' 



,(0)\AT,(0) 



(0)\ 



(*i°)|D|vl.(°))(vI/f |H| 



■E- 



EDM 



(0.6) 



E„ - Ei 



since both D and H^^^ are hermitian. In the above 
equation / represent intermediate states. 

It is obvious from the above equation that, the accu- 
racy of the calculation of Da depends on the excitation 
energies of the different intermediate states, the matrix 
elements of H|j!)'^ and D. 



approach, l^*!^^) implicitly contains all the intermediate 
states. 

The perturbed cluster operators can be written as 



Sy = 



(0.8) 



where T^^^ and Si, are the first order of corrections 
to the cluster operators T^^'' and si^^ respectively. The 
amplitudes of these operators are solved, keeping up to 
linear in EDM perturbed amplitudes, by the following 
equations 



(ci>S|<r(i) + H|-ps|$„) ^ 0, 
«|H!:'r(i) + H|^|$o) = 0, 



(0.9) 



and 



= -(<i>S|5W|a>.)ip, 



iKlK^Si'^ + (Hi^^T^i) + H|-PS ){i + 

= -«l4'^l'i>.>lP,(0.iO) 

where H'"^ is the Dirac-Coulomb (DC) Hamiltonian and 
H is defined as e~"^*°'He"^'''' which is computed after de- 
termining T^^') , IP is the ionization potential correspond- 
ing to the valence electron 'v' and the subscript N de- 
notes the normal form of an operator. We have used a, b.. 
and p,q.. etc. to represent holes and particles respec- 
tively, and are the single and double excited 
states respectively with respect to the \^v)- Using eqns. 
(0.3), (0.5), (0.8) and only keeping terms hnear in Gf, 
the expression for Da can be written as 



~ 1 +1^' _ 

< $.|^^^^'d(o)(i + 5^°^) + (1 + 5^°^')D(")5l^^ + g<°)'(r(i)^D(") + D(")r(i))5l°^ + (t^d^dw + D(")r(i))5^°^ |$„ > 



In the above expression we define D(o) = e^'°^ De"^'"' and N^^"^ = si^'' e"^*"' e^''°^ si'^'^ for the valence electron 
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TABLE L Excitation energies (cm El elements (a.u.) and magnetic dipole hyperfine structure constants (MHz) of important 
intermediate states in cesium and thallium. 





Cesium 


Thallium 


Initial state 
— > Final state 


6Sl/2 ^ 6pi/2 




6Sl/2 7pi/2 


6pi/2 7si/2 




6pi/2 8Si/2 


Excitation 

energy 

Experiment 


11229.38 
11177.84"^* 




21796.31 
21765. 30''''' 


26038.62 
26477.50° 




10462.32 
10520.0r 


IJjL LiaiiblLiOll 

amplitude 
Experiment 


4.53 
4.52(1)'= 




0.292 

0.284(2)'= 


1.84 
1.81(2)= 




0.57 


Atomic state 


6si/2 


6pi/2 


7pi/2 


6pi/2 


7Sl/2 


8Sl/2 


Hyperfine 
constant (A) 
Experiment 


2292.32 
2298.16'* 


284.86 
291.90(9)'' 


94.67 
94.35(4)'* 


21025.98 
2131T= 


11992.11 

12297"= 


4118.57 



References: " |21|1; '= (221; <* " 




(1) 



(1) " 



(c) 



^Iv 

(0) 



-(0) 



(d) 



'iv 



(1) 

Iv 



,(0) 



(g) 



(h) 



FIG. 1: Important Goldstone diagrams corresponding EDM 
amplitude calculations. 



'w' and each term is connected. 

The orbitals are constructed as linear combinations of 
Gaussian type orbitals (GTOs) of the form [l8l | 



(0.12) 



where k = 0,1,.. for s,p,.. type orbital symmetries re- 
spectively. For the exponents, we have used 



ai = ao/3* 



(0.13) 



We have considered ao as 0.00525 and (3 as 2.73 for all 
the symmetries for both the systems. To obtain Dirac- 
Fock (DF) wave function, we have undertaken 36si/2, 
34pi/2, 34p3/2, 324/2, 32^5/2, 3O/5/2, 3O/7/2, 20^7/2 
and 20(79/2 GTOs for cesium and 38si/2, 36pi/2, 36p3/2, 
36^3/2, 36^5/2, 30/5/2, 30/7/2, 25517/2 and 25519/2 GTOs 
for thallium respectively. All core (occupied) electrons 
are excited for both the systems. All orbitals are gener- 
ated on a grid IT^ using Fermi nuclear distribution |19| |. 



We have presented in table I properties related to the 
EDM in the ground states of atomic Cs and Tl. The 
square root of the products of the relevant hyperfine con- 
stants (A) provide an estimate of accuracy of the EDM 
matrix elements. The results of the various properties of 
Cs and Tl are in most cases in very good agreement with 
the measured values. 

TABLE II: Contributions from important RCC terms for 
Da/Cs calculations of 6s ■^S'1/2 and 6p ^Pi/2 states for ce- 
sium and thallium respectively in xlO^^^e — cm. 



Important 
terms 



Contributions from DF 

P,ttS-PS , ttS-PS pi 



Contributions from RCC 

Dr[i'+T(i)'D 

D(0)g<t' + S-'i'tDTO 
DWS^t'+s'i^+DW 

(1) , c(i)t?:?oTQ(o) 



j{0)t 

Norm, 



CV"Ht-)(0) c 



^2v ^ '-'2d 



Total 



Cesium 
6^ ^^1)^2 



Thallium 



6p2R 



(1) 

l/2_ 



-0.578 



3.217 



0.035 


3.056 


0.878 


4.453 


0.043 


-3.835 


0.015 


-0.304 


0.041 


0.174 


0.004 


0.023 


0.008 


-0.036 


0.019 


-0.032 


0.801 


4.056 



We have calculated the ratio of Da / Cs for Cs and Tl. 
Important Goldstone diagrams considered in these calcu- 
lations are shown in FIG. 1 and the contributions from 
significant RCC terms are given in table II. The lead- 
ing contributions for both Cs and Tl come from DS'^i^ 
(diagrams l((c)-|-(d))). These diagrams represent the 
DF, important pair correlation and a sub class of core- 
polarization effects. The DF contributions are -0.578 and 
3.217 in the units given in the table for Cs and Tl respec- 
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TABLE III: Comparison of Da/Cs results with others in 
xlO^^^e - cm. 





Cesium 


Thalhum 




Da/Cs 


Da/Cs 


Bouchiat [26] 


-0.689 




Venugopal [271 
Martensson-PendriU 


-0.805 




and Lindroth f28l 


-0.72(1±0.03) 


7±2 


Present 


-0.801(±0.004) 


4.056(±0.137) 



tively. Two different types of core-polarization effects 
represented by the terms HT^^ (diagrams l((a)-f(b))) 
and DS'i^ (diagrams l((e)-H(f)-H(g)-h(h))) make very 
large contributions but with opposite signs in the case of 
Tl. This is primarily due to the strong Coulomb and S-PS 
interactions between the Qpi/2 valence and the 6si/2 core 
electron. However, these interactions are much weaker in 
the case of the 6si/2 valence and the core electrons for 
Cs. The corresponding core-polarization contribution is 
therefore much smaller in size for Cs. It is of interest to 
note that the contributions of two other correlation ef- 
fects represented by the terms s'(l^'^DS'~^J and S'^^^^DS'^J,' 
are non-negligible. 

The results of various calculations of Da/Cs are pre- 
sented in table III. We have estimated the error by taking 
the difference between our RCC calculations with single, 
double and the leading triple excitations and just single 
and double excitations. Bouchiat 's calculation is based 
on an one electron potential and relativistic corrections 
are added to the wave functions Venugopal's result 
for Cs is obtained using a hybrid method which combines 
certain important features of many-body perturbation 



theory and the multi-configurational Dirac-Fock method 
[27l |. Martensson-PendriU and Lindroth have computed 
the Dirac-Fock and various types of core-polarization 
effects 01 . However, they have not considered pair- 
correlation effects which are included to all orders of the 
Coulomb interaction and one order in S-PS interaction 
in our work. 

We have obtained a new limit for the S-PS coupling 
constant Cs = (0.995 ± 1.756 ± 0.033) x 10"^ by com- 
bining our RCC calculation of Da/Cs for Tl and the 
measured value of Da for that atom |9|; the first uncer- 
tainty comes from the experiment and the second from 
theory. This is a significant improvement over the cur- 
rent limit Cs (2 ± 7) X 10"'^ |2||. It would be possible 
to improve this limit even further if the new generation 
of EDM experiments using cold cesium atoms reach their 
expected levels of accuracies 0, ■ The results of 

these experiments could then be combined with our cal- 
culation of Da/Cs for cesium (0.5% accuracy) to yield 
this new limit. 

In summary, we have calculated Da/Cs for Cs and Tl 
using RCC theory with accuracies of 0.5% and 3.3% re- 
spectively. Many-body effects were found to be of crucial 
importance, particularly for Tl. Our calculated value of 
Da/Cs for Tl in combination with the experimental re- 
sult of Da for the same atom gives the most accurate 
limit for the S-PS couphng constant to date. 
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ful communications. We would like to thank Professors 
Daniel Heinzen and David Weiss for information about 
the status of their experiments. B. K. Sahoo is thank- 
ful to DAAD for providing him with a scholarship. The 
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